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Abstract 


Although mutations drive the evolutionary process, the rates at which the mutations occur are 
themselves subject to evolutionary forces. Our purpose here is to understand the role of selection 
and random genetic drift in the evolution of mutation rates, and we address this question in asexual 
populations at mutation-selection equilibrium neglecting selective sweeps. Using a multitype 
branching process, we calculate the hxation probability of a rare nonmutator in a large asexual 
population of mutators, and find that a nonmutator is more likely to fix when the deleterious 
mutation rate of the mutator population is high. Compensatory mutations in the mutator 
population are found to decrease the hxation probability of a nonmutator when the selection 
coefficient is large. But, surprisingly, the hxation probability changes nonmonotonically with 
increasing compensatory mutation rate when the selection is mild. Using these results for the 
hxation probability and a drift-barrier argument, we hnd a novel relationship between the mutation 
rates and the population size. We also discuss the time to hx the nonmutator in an adapted 
population of asexual mutators, and compare our results with experiments. 


KEY WORDS: mutation rates, branching process, hxation probability, hxation time 
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Because most mutations are deleterious, the mutation rate can not be too high; in fact, in an 
infinitely large population, for a broad class of fitness functions, an error threshold has been s 


hown 


to exist above which t 

le deleterious effects of mutation 

1971; 

Jain and Krug 

2007 

). The mutation rate is not 


20071 1 ■ and it has 


been argued that the stochastic fluctuations in a finite population limit the evolution of mutation 
rates below a certain level since in small enough populations, the adv antage gained by lowering the 


2 OIOI ). Empirical data 


mutation rate can not compensate the effect of random genetic drift ()LynchI I! 
for organisms with widely different effective popul ation size shows a negative correlation between the 


deleterious mutation rate and the population size (jSuNG et al 


2013), and some quantitati ve insig ht 


into this relationship has been obtained by treating all deleterious mutations to be lethal (|Lynch 


20111 ). However, this is clearly an extreme scenario, and it is important to ask how the deleterious 


mutation rate evolves when mutations are only weakly deleterious. 

Many theoretical and experimental investigations have also shown that in an adapting asexual 


population, a mutator allele causing a higher mutation r ate t 


fixed (see a recent review by 


ran that of the nonmutator can get 


Raynes and SniegowskiI (1201411 1. Since the mutators produce not 


only deleterious but also beneficial mutations at a higher rate than the nonmutators, the mutato r 


allele can hitchhike to fixation with favorable mutations (Smith and Haigh 


19741 : 


Taddei et al. 


19971 1. However, once the population has reached a high fitness level, high mutation rates are 


is expected to decrease ( 

Liberman and Feldman 

1986 

). 

ndeec 

in some experiments 

(Trobner and Piecho( 

:;k] 

1984; 

Notley-Mc 

Tobb et al. 

2002; 

McDonald et al. 

2012, 


Turrientes et al. 


2013 


WiELGOSS et al 


2 OI 3 I ). the mutation rate of an adapted population 


carrying a mutator allele has been seen to decrease and the time to fixation has been measured, but 
a theoretical understanding of this time scale is missing. 

To address the issues discussed above, we study the fate of a rare nonmutator in a lar ge asexual 


population of mutators using a multitype branching process (|Patwa and Wahl 


important difference between the previous works on mutator hitchhiking (|Taddei et al 


2 OO 8 II. An 


Andre and Godelle 


200(tI : 


Wylie et al 


20091 : 


Desai and Fisher 


19971 : 


20 111 ! and our study is that here 


the mutator population is assumed to be at mutation-selection equilibrium and is therefore not 
under positive selection. However, compensatory mutations that alleviate the effect of deleterious 
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mutations are included in our model. We find that when only deleterious mutations are present, a 
nonmutator can get fixed with a probability that increases with the deleterious mutation rate of the 
mutator. Compensatory mutations in the mutator population are expected to decrease the fixation 
probability of the nonmutator, and we find that this intuition is indeed correct when deleterious 
mutations in the mutator are effectively lethal. But, surprisingly, when the deleterious mutations are 
mildly harmful, the fixation probability is found to initially increase and then decrease as the rate of 
compensatory mutations increases. Our study thus identifies the conditions under which the spread 
of nonmutators is suppressed in the absence of positive selection, and complements ea rlier works in 


which a mutator hitch 

likes with beneficia 

mutations to fixation ( 

Taddei et al. 

1997 

Andre and Godelle 

2006; 

Wylie et al. 

2009; 

Desai and Fisher 

2011 

)• 


Using our results for the fixation probability and a drift-barrier argument which states that the 
advantage offered b y a decrease i n the deleterious mutation rate is limited by random genetic drift in 


a finite population (jLYNCH 


20101 b we find that the deleterious m utation rate decre ases with 


increasing population size in accordance with experimental data (ISuNG et al 


2 OI 2 I ). H owever, unlf 


ce 


previous theoretical work that treats the deleterious mutations to be effectively lethal (jLYNCH 


2011 


he re we consider both strongly and weakly deleterious mutations, and not only reproduce the result 


m 


LynchI (|201ll b but also find a new scaling law in the latter case. We also use the results for the 


fixation probability to find the time to lower the mutation rate in an ad apted population o f 


mutators, and compar e our theoretical results with recent experiments (jMcDONALD et al. 


WiELGOSS et al. 


201a). 


2013; 


Model and methods 


We consider an asexual population in which the fitness of an individual with k deleterious mutations 
is given by W{k) = (1 — s)^, where the selection coefficient 0 < s < 1. A deleterious mutation is 
allowed to occur at a rate Ud and a beneficial one at a rate Ub < Ud- We are interested in the fate of 
a nonmutator that arises in this population and whose total mutation rate is smaller than that of the 
mutator. In a sufficiently large population of mutators in which stochastic fluc tuations due to 


genet ic drift may be ignored, this can be addressed using a branching process (jPATWA and Wahl 


20081 ). as described below. 
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The fixation probability TT{k,t) of a single copy of a nonmutator allele with fitness W{k) present at 


generation t changes according to (jJOHNSON and Barton 


2003) 


1 — 7r(fc, t) = exp 


^ k')7r{k',t + l) 


( 1 ) 


where W{t) = average fitness of the mutator population andp(A:,t) is the 

mutator frequency. The above equation expresses the fact that a single copy of the rare allele in the 
fitness class k whose offspring distribution is Poisson with mean W(k)/W(t) will be lost eventually if 
each of its offspring, which may undergo mutations with probability M{k ^ k'), do not survive. 

’ iere we consider strong mutators whose muta tion rate is much higher than that of the nonmutator 


SNiEGOVySKi et al 


19971 : 


Oliver et al. 


2000 ), and therefore neglect the mutation rate of the latter 


in most of the following discussion (however, see Fig. [T|). We also assume that the mutator 
population is at mutation-selection equilibrium as is likely to be the case in large populations that 
have been evolving for a long time in a constant environment. As a result, the probability 7r{k, t) 
becomes time-independent. These considerations lead to a relatively simpler, but still highly 
nonlinear equation given by 


1 — TT{k) = exp 


W{k) 7r(fc) 

w 


( 2 ) 


The above expr ession, of course, reduces to the well known single-locus equation (jFiSHER 


Haldane 


1922 


19271 1 when the nonmutator can be present in only one genetic background, but here we 


are dealing with a multitype branching process since a nonmutator can arise in any fitness class. 
he total fixation probability is obtained on summing over all genetic backgrounds 


(Johnson and Barton 


2003), 


VTtot = ^p(fc) 7r(/c) 


(3) 


k=0 


where the probability that a nonmutator arises in a background of k deleterious mutations is given 
by the mutator frequency p{k) in the stationary state. 

Although the steady state frequency p{k) in the abs ence of compensatory mutations that mitigate 


the effect of deleterious mutations is known exactly (IKimur A and Martjyama 


19661; 


Haigh 


1978|), 


the corresponding solution with nonzero Ub is not known. We therefore compute the mutator 
frequency numerically for nonzero Ub using (jA.ljl given in Appendix [Aj and use these results in ([2]) 
to find the fixation probability for arbitrary Ub- To make analytical progress, we use a perturbation 
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theory in which the effect of the small dimensionless parameter Uh/s can be studied by expanding 
the quantities of interest in a power series in C/ft/s, and write 

n=0 ' n=0 ^ 

The terms '7ro{k) and poik) corresponding to n = 0 in the above expansion give the results in the 
absence of compensatory mutations, and in Appendix we calculate the stationary state fraction 
p{k) to linear order in Uh/s. 


Results 

FIXATION PROBABILITY 


In the absence of compensatory mutations: We first consider the case when Uh = 0. Taking 
the logarithm on both sides of dS]), and expanding the left hand side (LHS) up to vrg(/c), we find that 
either 'iro{k) = 0, or 


TToik) = 2 


V 


2s{ko - k) , 


(5) 


where the average fitness V 


anc 

the 

( Kimura and Maruyama 

1966; 

Haigh 

1978) 


an d the average number of deleterious mutations /cq = U^/s 


1973)- The last expression on the right hand side (RHS) of 


([5]) is obtained by expanding the exponentials as the parameters and s are small. Since the 
fixation probability must not be negative, the expression ([5|) is valid when k < [fcoj, and the solution 
7 ro(fe) = 0 holds otherwise. Here [xj denotes the largest integer less than or equal to x. More 
generally, a nonmutator can get fixed if its fitness W(k) ~ is larger than the average fitness 
e-sk q£ mutator populatio n, or k < [k \, k being the average number of deleterious mutations 


(Johnson and Barton 


2002 ). 


Equation ([5|) shows that the fixation probability vro(fe) decreases as the number of deleterious 
mutations increase, as one would intuitively expect. However, the probability po{k) that a 
nonmutator would arise in a background with k < ko deleterious mutations increases. On summing 
over the backgrounds in which a nonmutator can arise, as explained in Appendix [HI we find that the 
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total fixation probability falls in two distinct regimes defined by whether Ud is below or above s: 


[ko] 

TTo = '^'n:o{k)po{k) 
k=0 


= 


2Ud , Ud<^s 

V^, u^»s 


( 6 ) 


For fco <C 1, since a mutation is costly, it can be treated as effectively lethal ([Johnson 


199911 . In this 


situation, the advantage conferred by the nonmutator is simply given by 1 — e ~ Ud and the 


classical result f or the single locus problem gives the fixation probability to be 2Ud ([Fisher 


Haldane 


1922 : 


19271 ). For A:o S> 1, the total hxation probability apparently receives contribution from ko 


genetic backgrounds, but merely \/ko genetic backgrounds are actually relevant because the 
Poisson-distributed frequency po{k) has a substantial weight for fitness classes that lie within a 
width y/ko of the mean (also, see Appendix [Bj) . Equation ([6|) shows that for fixed s, the nonmutator 
is more likely to be hxed when Ud is large. But, for a given Ud, the hxation probability initially 
increases with the selection coefficient and then saturates to 2Ud- In Fig. [H the analytical results 
above are compared with those obtained by numerically iterating ([2|) and ([T|) when the mutation rate 
of the nonmutator is zero and Ud/50 respectively, and we see a good agreement in both cases. 
Including compensatory mutations: We now study how compensatory mutations in the mutator 
population affect the hxation probability of the nonmutator. Figure [2] shows that when ko <C 1, the 
hxation probability decreases with Uh, but for fco ^ 1) h changes nonmonotonically: it hrst increases 
and then decreases with increasing Ub- To understand this behavior, consider the change 

= TTtot — ttq in the hxation probability due to compensatory mutations which is simply given by 


Lfcj 

(^vTtot = '^Po{k)S7r{k) + TTo{k)5p{k) -h 5p{k)5'K{k) . (7) 

k=0 

When Ub is nonzero, the change in the hxation probability 5Ti{k) = 7r(fc) — TTo{k) and the mutator 
frequency 6p{k) = p{k) — Po{k) behave in a qualitatively different manner. With increasing Ub, the 
average fitness of the mutator population increases which, by virtue of ([2|), decreases the hxation 
probability of the nonmutator, i.e. J7r(fc) < 0. However, since the frequency of individuals with less 
deleterious mutations increases when Ub is nonzero, the change in the mutator fraction 5p{k) > 0. 
Thus the change in the total hxation probability given by ([7|) receives both positive and negative 
contributions, and it is not obvious which one of these factors would have a larger effect. 

To address this question, we calculate the hxation probability for small Ub/s as described below. 
Substituting ([1|) in the expression ([7|) for Jvrtotj and neglecting terms of order {Ub/s)"^ and higher, we 
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find that ^vrtot ~ {Ub/s)T^i, where 


\ko\ 

TTi = '^Po{k)7ri{k) +pi{k)7ro{k) . (8) 

k=0 

The contribution Tri{k) is calculated in Appendix [Cl and we find that 

7ri{k) ^ -2sko{l - Tro{k)) , k < [ko\ , (9) 

which is negative, as expected. An expression for the fraction pi{k) is obtained in Appendix |Al and 
its behavior is shown in Fig. [3] for small and large fco- For small ko, the frequency po{k) is close to 
one in the zeroth fitness class and zero elsewhere. But the correction pi{k) is negligible in all the 
fitness classes. For large fcoi the contribution pi{k) is significantly different from zero in many fitness 
classes and can be approximated by 

pi{k) = koPo{k) In , A: > 1 . (10) 

Thus, as claimed above, the fraction pi{k) is positive for k < ko and negative for k > ko (also, see 
Fig. [3]). 

When Ud <C s, as already mentioned, the fraction pi{k) is negligible in all the fitness classes and 
po(0) « 1. Using these results in (l8|) and (|9]), we get vri = —2sko, and thus 

— ) Ufj < Ud < s. (f-f) 

TTo S 

This reduction in the fixation probability of the nonmutator when Uf, is nonzero is expected since the 
effect of compensatory mutation is to restore the mutators that have suffered lethal mutation to the 
zeroth mutation class thus enabling them to offer competition to the nonmutators. 

When Ud S> s, as shown in Appendix [Cl we can obtain a quantitative estimate of the initial increase 
in (fvTtot by calculating the sum on the RHS of ([8|) to obtain (|C.5p . and thence 

7i > < S < Ud ■ (f-2) 

TTo 2s 

Thus we find that for small Ub, the increase of the mutation frequency in fitness classes with fewer 
deleterious mutations dominates the increase in the mutator fitness resulting in positive (Jvrtot- 
However, for large Ub, the net change in the fixation probability is negative since the last term in the 
summand of ([7|), which is also negative, enters the picture. Since the maximum in (^TTtot occurs at 






large Uh/s, the perturbation theory described here can not capture the eventual decrease in this 
parameter regime. A quantitative comparison of the results obtained by numerically solving ([2]) and 
(lA.ll) for arbitrary with the analytical results (1111) and (1121) for small C/fe/s is shown in Fig. [21 and 
we observe a good match when Uh/s is small. For large Uh/s and Ud/s, a fit to the numerical data 
shows that the fixation probability decreases linearly with Ub- 


EVOLUTION OF MUTATION RATES IN FINITE POPULATIONS 

The drift-barrier hypothesis states that in a finite population, the beneficial effect of lower 
deleterious mutation rate can be outw eighed by the stochastic effects of random genetic drift which 


limits the evolution of mutation rates (Lynch 


20101 1. In a finite population of size N, a mutation 


that decreases the deleterious mutation rate confers an indirect selective advantage and will spread 
through the population. However, as Ud decreases, the fixation probability of such a mutant 
decreases until it reaches its neutral value vTneu = 1/N. Here we have calculated the fixation 
probability vro neglecting stochastic fluctuations. The full fixation probability H that includes the 
neutral and the large population limit may be obtained as follows. 


The fixation time for a mutator in a finite population of n onmutators when a . 
deleterious has been calculated using a diffusion theory bv i.lAiN and NagarI i 


1 mutations are 


201j), and shown to 


increase exponentially with the population size. The fixation probability ^ e is thus 


exponentially small in the population size (jKlMURA 


1980 


Assaf and MOBILIA 


20111 ). where we 


have identified the rate of decrease of fixation probability with a selection coefficient 25. This 
effective selection coefficient is found to match exactly with the result Q for the fixation probability 
ttq obtained here using a branching process. Although this is not a rigorous proof, these observations 
strongly suggest t hat the fixation probability of a nonmutator in a finite population of size N is of 


the classical form (jKiMURA 


1962) 


H = 


1 — e 


-25 


1 - e-2^5 ’ 

where S = 7ro/2. We also mention that the probability 25 depends on the difference in the 


(13) 


deleterious mutation rate of the mutator and t 


re nonmutator when the mutation rate of the 


nonmutator is nonzero ( 

Jain and Nagar 

2013 

), and has also been shown to be insensitive to the 

distribution of selective effects ( 

Desai and Fisher 

2011 

)• 
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Thus, according to (1131) . a crossover between positive selection and neutral regime occurs when 
TTo ~ N~^ and gives a lower bound on the mutation rates. We recall that the fixation probability tto 
in ([S]) shows a transition when ~ s, and at this mutation rate, the fixation probability vro ~ s. 
This translates into a change in the behavior of Ud when Ns crosses one, and we have 


Ud ~ < 


N-^ , > 1 . 


(14) 


Thus in the weak selection regime {Ns <C 1), the deleterious mutation rate depends on the selection 
coefficient, and decreases faster than when the selection is strong. Figure U] shows the preliminary 
results of our numerical simulations for a finite size population of mutators with mutation rate Ud in 
which nonmutators with mutation rate Ud/2 can arise with a certain probability. This population of 


nonmutators and mutators evolves via stanc 

ard Wright-Fisher dvnamics, and the time to fix the 

nonmutators is measured ( 

Jain and Nagar 

2013 

). For a fixed N, the fixation time is found to 


increase as the mutation rate of the mutator is decreased until a minimum mutation rate is reached 
below which the hxation time remains constant. This lower bound, shown in Fig. SI exhibits 
different scaling behavior in the weak and strong selection regimes, in accordance with (I14|) . 


Discussion 

Fixation probability: A rare mutator arising in a population of nonmutators carries a higher load 
of deleterious mutations but offers indirect benefit by producing more beneficial mutations. The 
fixation probability of a rare mut ator in a finite nonm utator population has been studied by 


Andre and Godellh (j2006ll and 


Wylie et al. 


(j2009|) analytically, and found to vary 


nonmonotonically with the mutation rate of the mutator. It has been shown that the fixation 
probability is of the classical form ()13p where the effective selection coefficient S when scaled by the 
selective advantage s increases (decreases) when the ratio of mutation rate to selection coefficient is 
below (above) one. Here, we studied a situation in which a nonmutator appears in a mutator 
population and is benehcial since it produces fewer deleterious mutations, and calculated its fixation 
probability vrtot using a branching process. The mutator population is assumed to be at 
mutation-selection balance and therefore, by definition, selective sweeps resulting in the spread of 
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favorable mutations are neglected. However, it is interesting to note that the scaled fixation 
probability of the nonmutator obtained here also changes its behavior when the deleterious mutation 
rate is o f the order of the selection coefficient, see ([S]). Our work significantly extends the previous 


result of 


Lynch! ([20111) as the deleterious effect of mutations is allowed to be mild here, and 


therefore we are dealing with a truly multilocus problem. 

Compensatory mutations that alleviate the effect of deleterious mutations are found to have a 
surprising effect on the fixation probability of the nonmutator. Although they improve the fitness of 
the mutator population, it also means that the nonmutator can arise in a better genetic background 
where it has a better chance of fixation. Thus, compensatory mutations affect both the resident 
mutator population and the invading nonmutator allele in a positive manner. The effect of these two 
factors on the fixation probability of the nonmutator is, however, opposite and can result in an 
unexpected increase in the fixation probability of the nonmutator when compensatory mutations are 
present. Here we have shown analytically that this scenario is realised when the mutations are 
weakly deleterious and the compensatory mutation rate is small, as illustrated in Fig. |2l The 
increase in the fixation probability due to compensatory mutations can be quite high, but we do not 
have analytical estimates for this. An exact solution of (jA.ip would, of course, pave the way for a 
better analytical understanding but is currently not available. 

Fixation time: In a maladapte d asexual population, the mutat ors can sweep the population since 


they facilitate rapid adaptation ( Raynes and Sniegowski 


2014lf . But as the population adapts and 


the supply of beneficial mutations diminishes, mutators have a detrimental effect on the population 
fitness and a mutation that lowers the mutation rate is favored. In bacteria E. coU, several genes 
(such as mutT and mutY) are involved in avoiding or repairing the errors that oc cur duririg the 


replication process, and defects in these genes can lead to the mutator phenotype ([Miller 


199®. 


But compensatory mutations in the defec tive error-repair machinery can reduce the mutation rate, 


at least, partially ( WiELGOSS et al 


2013]). We therefore model this situation by assigning a 


probability b with which mutators can convert into nonmutators due to a mutation in the 
proofreading or error-repair region. In E.coli, the conversion probability f from nonmutat or to 


mutators has been estimated to be ~ 10 ® per bacterium per generation ( Boe et al. 


2000). But the 


probability b for the reverse m utation is not known, a lthough one expects b < f, possibly because it 


is a gain-of-function mutation ([Wielgoss et al. 


2010 ). 
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When the rate Nb at which the nonmutators are produced from the mutators is smaii enough that 


the new aiieies behave independentiy, the time taken to fix the nonmutator popu 


T = (NbTTiot) In a iong term evoiution experiment on E.coli, 


WiELGOSS et al. 


at ion is given by 


( 2013 ) found the 


mutation rate to decrease by about a factor two in a neariy adapted mutator popuiation with a 
mutation rate 150 times that of the wiidtype in t wo lineages. As th e population size in Lenski’s 


experiments has been estimated to be about 10^ (|Wahl et al. 


2 OO 2 I'). the product Nb can be at most 


ten which is not too large. We first note that in the experiment of 


WiELGOSS et al. 


(I 2 OI 3 ), the 


fixation time was longer in the lineage in which the mutation rate decreased by a smaller amount, in 
accordance with ([6]). To make a quantitative comparison, we consider the ratio of the times for the 


two lineages, as T dep ends strongly on t 


the data in Table 2 of 


WiELGOSS et al. 


re pr obability b which is not known experimentally. Using 


(120131 1. we find the ratio of fixation time in mutT mutY-L 


background to that in mutT mutT-T'background to be 9209/5157 « 1.8. The theoretical formula dH), 
on replacing by the difference between the mutation rate of the nonmutator and mutator yields 
1.5 (1.2) when mutations are assumed to be strongly (weakly) deleterious and the selection coefficient 
same in both lineages. Since ([6]) is obtained assuming that the mutators are strong whereas the 
mutation rates decreased merely by a factor two in the experiment, a more careful examination is 
needed. Solving ([T|) numerically in the stationary state, we find that the ratio is unaffected whe n the 


mutat ions are strongly deleterious. But using the mutation rates in Table 2 of 


WiELGOSS et al. 


(j2013ll and s ~ 0.01 yield the ratio to be about 4.5. Although the theoretical conclusions (1.5 — 4.5) 


are in reasonable agreement with experiments, the above analysis suggests that the reversion 
probability b may not be too small {i.e. N b >, 1), and a more sophistic ated theory that takes care of 


the interference between the nonmutators ( Gerrish and Lenski 


19981) may be required to obtain a 


closer match. We close this discussion by noting that in an experiment on S. cerevisiae in which the 


adapted population reduced its genome-wic 


experimental lines (IMcDonald et al 


e mutation rate by almost a factor four in two of the 


201211 . the fixation time seems to increase with the mutation 


rate, in contradiction with the experiment of 


WiELGOSS et al. 


((20131) and the theory presented here. 
Evolution of mutation rates: Exp eriments show tha t the mutation rate decays as for 


prokaryotes and N for eukaryotes (|SUNG et al. 


2 OI 2 I I. The population size and deleterious 


mutation rates are negatively correlated since deleterious mutations can get fixed in small 
populations due to stochastic fluctuations, but not in large populations where the genetic drift is 
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ineffective ( Lynch 


2010|). Here, we have shown that a reciprocal relationship between the population 


size and mutation rate holds for large populations, but for small populations, the deleterious 
mutation rate decreases much faster, see Fig. [H This is in contrast to experimental results 
mentioned above where the data has been htted assuming a single scaling law. In view of our 
theoretical results discussed above, a more careful analysis of experimental data is required. 

While the evolution of deleterious mutation rate has received much attention, to the best of our 
knowledge, analogous theoretical predictions for the benehcial mutation rate are not available. As 


large populations experience clonal interference (jGERRiSH and Lenski 


19981 1 which results in the 


wastage of benehcial mutations, the r ate of benehcial mutati ons is observed to be smaller in large 


populations in microbial experiments (|Perfeito et al. 


20071 1. An understanding of the relationship 


between the population size and the rate of benehcial mutations would be an interesting avenue to 
explore. Other potential factors that can affect the correlation between the mutation rate and the 
population size include epistasis and recombination. Here, we have also ignored the c ost of hdelity, 


and it remains to be seen how the r esults presented here are affected on including it (IKimura 


Kondrashov 

1995; 

Dawson 

19981 


1963; 


19981 1. A more detailed understanding of the mutation rates, both 


empirically and theoretically, remains a goal for the future. 
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A Mutator frequency when compensatory mutations are included 

For small selection coefficient and mutation rates, the mutator frequency p{k,t) obeys the following 
continuous time equations: 

dp{0, t) 


dt 

dp{k,t) 

dt 


= -Udp{0, t) + Ubp{l, t) + skp{0, t) 

= -Up{k, t) + Udp{k - l,t) + Ubp{k + l,t) - s{k - k{t))p{k, t) , 


(A.l) 


where U = Ud + Ub and k{t) = YlT=o ^ stationary state, the LHS is zero and the 

frequencies are time-independent. On dividing both sides of the above equations by s, we find that 
the stationary frequency p{k) depends on the ratios Ub/s and Ud/s. We first expand the fraction 
p{k) in a power series about [4/s = 0 as 

'Ub'''^ 


P{k) = 


n=0 


Pn{k) , 


(A.2) 


where Pn{k) is proportional to the nth derivative of p{k) with respect to Ub/s evaluated at Ub = 0. 
The lowest order term po{k) is the solution of the steady state of (lA.ip in the absence of 


comp ensatory mutations, and is known to be a Poisson distribution with mean ko = Ud/s (IHaigh 

lOTsI): 


Po{k) = , k = 0,l, 


(A.3) 


To find the solution with nonzero Ub, we first set the LHS of (jA.ip equal to zero, and substitute 
(]A.2I) in these equations. On neglecting the terms of order (Ub/s)"^ and higher, we obtain the 
following equations for pi{k): 


kiPo{0) = -po(l) (A.4) 

kopi{k - 1) - kpi{k) + hpoik) = po{k) - poik + 1) , k = 1,2, ... (A.5) 

where ki = YlT=okpiik). Equation (IA.4D above immediately yields ki = —k^. Thus, as expected, the 
effect of compensatory mutations is to decrease the deleterious mutations in a population. Using this 
result in (jA.Sp . after some simple algebra, we get the following one-term recursion equation for 
Pi{k) , k>l: 

Pi{k) = ^ Pi{k “ 1) “ Q + 'Pdk) , (A.6) 

which can be iterated easily to give 

Pi{k) =-^ Pi{U) - Po{k) [ko{Hk+i - 1) + Hk] , (A.7) 


14 

















where the harmonic number the fraction pi(0) is determined using the 

normalisation condition, viz., YlT=oPi^') ~ Since the fraction po{k) already satisfies this condition, 
we have the constraint YlT=oP^i^) ~ using which, pi(0) can be found. For large k, using 
Hk « In/c in (|A.7|) . we obtain the expression (fTOl) . 


B Fixation probability in the absence of compensatory mutations 

To find the total hxation probability given by ([3]) , we use the expression ([5|) for the fixation 
probability and (IA.3I) for the mutator fraction po{k) which is a Poisson distribution with mean ko. 
When ko <C 1, we have vro ~ po(0)'/ro(0) = 2Ud- But for ko ^ 1, on summing over the mutator 
backgrounds in which a nonmutator can arise, we obtain the total fixation probability to be 




[feoj 

'^Po{k)TTo{k) = 2s 
k=0 


e-L^oJ L^ojlfcoJ+i 


(B.l) 


On using the Stirling’s formula xl ~ y/27rx{xjey for large x in the last expression, we immediately 
obtain Q. Another way of seeing the result in the large ko regime is by approximating the 
Poisson-distributed po{k) by a Gaussian with mean and variance equal to ko and thus obtain 

fko _ I (k-kgf rl 2 

TTo ~ 2s / dk{ko — k) e ^*0 ~ 2s\Jko / dx xe ^ (B.2) 

y ko Jo 

where we have used the fact that the mutator frequency is substantial in the fitness classes lying 
within a distance of the mean. 


C Fixation probability when compensatory mutations are included 

Inserting TTfot = vro + {Ub/s)TTi and W = Wo + {Ub/s)Wi in ([2]), and using the exact equation for 
TTo (A;), we get a rather involved expression for 7ri(A:) given by 

Tf/O- 

Wo Wo-W(A:)(l-7ro(A:)) ’ 

Since all the parameters are smaller than one, we work with the approximate expression ([5|) for the 
probability 7ro(fc) and arrive at ([9]). 
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We now calculate the contribution tti given by ([8]) when fco !:$> 1 using the expression ([9]) for 'ni{k) 
and the frequency pi{k) in (fT0|) . We have 

C» CX) OO 

TTi = —2sfco(l ~ E PQ{k)) + Is'^iko - k)pi{k) - ^ pi{k)'KQ{k) (C.2) 

fc=[fcoJ A:=LfcoJ 

OO OO 

= 2skQ ^ po{k)-2 ^ {ko - k) spi{k) , (C.3) 

fc=[fcoJ fc=[fcoJ 


where we have assumed that Ub, Ud, s are small, but Ub/s and Ud/s are finite. The last expression is 
obtained on using the normalisation condition '^’^^oPiik) = 0 and the expression for the average ki. 
For large ko, we approximate the Poisson distribution poik) by a Gaussian as 

1 / h — hr. \ 

Approximating the sums in (IC.Sp by integrals, we finally have 


'Hi 

s 



2ko 



(C.5) 


where we have carried out an integration by parts in the second integral on the RHS and neglected 
subleading terms in ko- 
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Figure 1: Dependence of the fixation probability obtained using a multitype branching process on the 
deleterious mutation rate for two values of the selection coefficient s and compensatory mutation 
rate Ub = 0. The points are obtained by numerically solving ([2]) when the mutation rate of the 
nonmutator is zero (O, □), and the stationary state solution of ([I]) when the nonmutator’s mutation 
rate is 50 times lower than that of the mutator (+, x). The lines show the analytical result (l6|). 
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Figure 2: Dependence of the fixation probability obtained using a multitype branching process on the 
compensatory mutation rate Uh for two values of Ud/s. The points show the numerical solution of ([2]) 
and the lines show the analytical results m and (jl2l) . The broken curve for Uh/s > 0.1 is a linear 
fit, 0.1 — 0.24:Ub/s, to the numerical data. For Ud/s = 0.1, the ratio Ub/s is also below 0.1 since Ub is 
assnmed to be smaller than Ud- 
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Figure 3: Change in the mutator frequency when compensatory mutations are included, 6p{k) = 
p{k) — Po{k) for Ub = 10“^. The points (O, □) are obtained by numerically iterating (lA.ip and (+, x) 
show the perturbation theory result (jA.7jl . and we observe a good agreement. The simple expression 
(fTOl) for large Ud/s is also shown (lines). 
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Figure 4: Relationship between the deleterious mutation rate and the population size for selection 
coefficient s = 0.0064 when compensatory mutations are absent. The points are obtained by numerical 
simulations of a Wright-Fisher process, and the lines show the iV-dependence in (jl4p . 
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